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Abstract We use a simple model for biological ageing to study the mor- 
tality of the population, obtaining a good agreement with the Gompertz law. 
We also simulate the same model on a square lattice, considering different 
strategies of parental care. The results are in agreement with those obtained 
earlier with the more complicated Penna model for biological ageing. Finally, 
we present the sexual version of this simple model. 
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1 Introduction 

One theory for biological ageing is the accumulation of bad inherited muta- 
tions. While the well known Penna model [IJ uses a genome represented by a 
bitstring, this simple model reduces the genetic information, transmitted 
with mutations from one generation to the next, to two integers, the mini- 
mum reproduction age a m and the genetic death age aa- At each iteration, 
every individual with age a between a m and ad gets at most one child (birth 
rate 6 = 1) which inherits these numbers, apart from random changes (mu- 
tations) by ±1. Individuals die with certainty if their age reaches a d , and 
they die earlier with probability N(t)/N max at every iteration, where N(t) 
is the actual population and N max a so-called carrying capacity taking into 
account the limitations of food and space (Verhulst factor). The probability 
Pb(i) for individual i to get a child takes into account the well known (see 
e.g. 0) trade-off between fecundity and survival and is the lower the larger 
the difference a d — a m is: 

p b (i) = (1 + 0M)/(a d ii) - a m (i) + 0.08) 



1 



where the constant 0.08 ensures that the denominator is never zero. (For the 
mutations we require < a m < < 32 to facilitate comparison with the 
Penna model which typically has 32 age intervals.) 

This simple model allows for a self-organization ( "emergence" ) of a broad 
distribution of a m values, and a relatively more narrow distribution of ad 
values (similar but not identical to those shown below in Fig. 6; the a m 
distribution gets more narrow relatively if, as in section 4, an increase of a m 
is coupled to an increase of a^.). It explained the catastrophic senescence of 
Pacific Salmon ||, the vanishing of cod in the northwest Atlantic through 
overfishing [0] , and allowed to take into account the social needs of humans 
for a minimum population size ||. But only with a birth rate increasing in 
old age U could a reasonable agreement with the Gompertz law be achieved. 

In the next section we present the mortality obtained with a slightly 
modified version of this model, which then agrees with the Gompertz law. In 
section 3 we present the results of the model when the individuals are placed 
on a square lattice and live under some parental care strategies. Finally, in 
section 4 we present the sexual version of the model and conclusions. 



2 The Gompertz law 

To get the Gompertz law of an exponentially increasing mortality function 
of age a, we now follow J7J and assume the genetic deaths to be probabilis- 
tic instead of deterministic, with probability pd = 1/{1 + exp[(3(ad — a)}} 
and (3 = 0.7; the deterministic case is recovered for (3 — > oo. The sur- 
vival probability thus decays like a Fermi function in physics. Fig.l shows a 
good agreement with the exponential increase of the mortality function || || 
/i(a) = \n[S(a — l/2)/5 , (a + l/2)] for humans at middle age (where S(a) is the 
number of survivors at age a; fi approximates the derivative —d(hiS)/da). 
For the deterministic case, (3 — > oo, the analogous plot is strongly curved. Al- 



ternatively one may omit [10] the condition that no age exceeds 32, increase 
the constant from 0.08 to about 2, and set (3 = 0.3; then Fig. 2 shows a similar 
increase of the mortality function with age, with a reasonable distribution of 
genetic death ages which is no longer cut off at 32. 
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3 Parental care on a square lattice 



Now we return to the original deterministic version of the model and place the 
individuals on a square lattice. The carrying capacity N max is now replaced 
by a maximum occupation per site N max (i,j), such that the new Verhulst 
probability for an individual to die due to environmental restrictions depends 
on the current number of individuals N^^(t) on each site. We give to each 
individual a probability p m to move to the neighbouring site that presents 
the smallest occupation, if this occupation is also smaller or equal to that of 
the current individual's site. The strategy of child-care consists in defining a 
maternal care period A pc during which no child can move. We considered the 
following conditions: (a) if the mother moves, she brings the young children 
with her; (b) the mother cannot move if she has any child still under maternal 
care. In Fig. 3 we show the lattice configurations after 800,000 steps for case 
(a). The result of case (b) is very similar and not shown. These same 
strategies were simulated before using the Penna model | TI| , giving these 



same results. Also following |TT| , we simulated a third case which corresponds 
to case (a) with the restriction that if the mother dies, her children still under 
maternal care have a further probability to die. Figure 4 shows the spatial 
configuration for this case (c). Again the result is the same obtained with the 
Penna model, despite of the simplicity of the present model. Fig. 5 compares 
the population sizes of the three cases. The parameters used to produce 
figures 3, 4 and 5 are: initial population = 10,000 individuals; maximum 
occupation per site = 30 individuals; initial minimum reproduction age a m = 
1; initial genetic death age a d = 16; birth rate 6 = 2; maternal care period 
A pc = 2 and lattice size = 150 x 150. 



4 The sexual version 

Now the population consists of males and females. At each iteration, every 
female with age a between a m and a d randomly chooses a male, also with 
age between a m and a d , to mate. The child randomly inherits the a m and 
ad values from one of the parents, independently. Mutations are introduced 
in the following way: two random numbers are generated; if they are both 
positive, a m and a d are both increased by one. If they are both negative, a m 
and a d are decreased by one; if the signs of the random numbers are different, 
no mutation occurs. Notice that in this way the difference a m — a d is main- 
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tained constant. This strategy simulates the fact that a recessive mutation, 
to be effective, must be inherited from both parents; alternatively, it can be 
interpreted as antagonistic pleiotropy, the trade-off between fecundity and 
longevity ||, |3|. One may also regard this sexual version as an alternative 
to the requirement of a minimal population for human society ||: a single 
individual cannot reproduce. In Fig. 6 we show the histograms of a m and a^, 
that is, the number of individuals with a given value of and the number 
of individuals with a given value of a m . It can be seen that there is also, 
as in the asexual case @, a self-organization of distributions of a m and a^. 
We have also tried to introduce the mutations in a m and ad independently 
of each other: Then no self-organization appears, that is, the a m value goes 
always to 1 and the value, to 32. 



As in many previous simulations, starting with Redfield [ 12], [0|, the 



sexual version does not give a clear advantage over the asexual one and 
additional assumptions like in |L4[ may be needed to justify sex. 



5 Conclusion 

The simple alternative 0] to the more widespread Penna model [|l], || of bi- 
ological ageing via inherited mutations was modified here to give reasonable 
agreement with the Gompertz mortality law at middle and old age, to allow 
a more realistic description of spatial fluctuations on a lattice, and to include 
sexual reproduction. In all cases the results were similar to those obtained 
earlier in the Penna model. Further research could concentrate on the sim- 
plification to omit the restriction ad < 32 and to omit the parameter 
now taken as 0.08. 
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-d[ln(survivors)]/da for 9.2 million individuals with (+) and without (x) Verhulst deaths 
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Figure 1: Mortality function // versus age (in iterations) giving rough agree- 
ment with the Gompertz law; (3 — 0.7 with and without Verhulst deaths 
using 9 million individuals. 
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600 million individuals, no maximum age, eps=2, beta=0.3, includes Verhulst deaths, t up to 20000 
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Figure 2: As Fig.l but (3 = 0.3 with 609 million individuals including 
Verhulst deaths, no maximum age and 0.08 replaced by 2.0 in the birth rate. 
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Figure 3: Final distribution of the population on the lattice considering 
that if the mother moves, she brings the children under maternal care (age 
< 2) with her (case a). The final distribution of the lattice considering that 
the mother cannot move, if she has children under maternal (case b), looks 
very similar (not shown). 
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Figure 4: Final distribution of the lattice considering that if the mother 
moves she brings the children under maternal care with her, and that if she 
dies, her children have a probability 0.9 to die in their first year of life, and 
0.3 in their second year of life (case c). In this case we used a maximum 
occupation per site equal to 42 individuals. Even so, the final configuration 
of the lattice is much smaller than those obtained in cases (a) and (b). 
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Figure 5: Number of individuals as a function of time for cases: a (circles), 
b (triangles) and c (squares). 
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Figure 6: Distribution of minimum age of reproduction (circles) and of 
genetic death age (triangles) after a stationary state has been reached for 
the sexual population. Parameters: initial population = 25,000 males plus 
25,000 females; maximum population size N max = 500, 000, initial a m — 1 
and initial ad = 16. 
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